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Motivations 


From “Deposition With Hot Streaks in an 
Uncooled Turbine Vane Passage’, B. Casaday, et 
al., J. Turbomach, 2013 Vol. 136 (Permission from 
Prof. Bons and thanks to Dr. Mike Dunn @ OSU) 


S 


Designing high-pressure turbines 
(HP Ts) for peak temperatures at 
the combustor exit =» More 
cooling air Less cycle efficiency 
Designing HPTs for the mean exit- 
temperature at the combustor exit 
= More local hot spots (hot 
streaks) = Less gas turbine 
durability 


CFD should give some design 
guidelines 


Important to understand core engine component 
interactions, such as combustor-turbine interactions 


Energy Efficient Engine— GE design, 80s - eS 
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High pressure turbine (HPT) 


E% is a double-annular and compact combustor, intensively investigated in 
the mid-1970s and 1980s and set a historically important milestone toward 
more fuel efficient jet engines (and now geometry is publically available!) 
HPT for the GE E8 is a two-stage, low thru-flow design for moderate loading. 
Relative position of the stator with respect to the fuel nozzle, so-called 
~ clocking’, is Critical. 


Features of Open National Combustion Code 
(OpenNCC) 


¢ OpenNCC is the releasable version of the National Combustion 
Code (NCC), which has been continuously updated for more than 
two decades at NASA Glenn Research Center (GRC) 


¢ Main Features 
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Numerics: Jameson-Schmidt-Turkel (JST) scheme and Roe's upwind 
scheme, and Advection Upstream Splitting Method (AUSM)(*-3.17) 


Turbulence: Cubic non-linear k-¢4) model with the wall function, Low- 
Re model 


Combustion: Finite Rate Chemistry, low dimensional manifold with 
EBU, PDF, Linear Eddy Model (LEM)‘) 

Spray: Lagrangian liquid phase model(®®) 

Other features: Low-Mach preconditioning, transition model), 
unstructured mesh, adaptive mesh refinement (AMR), massively 
parallel computing (with almost perfectly linear scalability achieved for 
non-spray cases up to 4000 central processing units) 
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Problem Setup 


v Geometry (clocking) 
v¥ Mesh 
Y Numerical Setting 


Y Boundary Condition and Operating Condition 


24 Degrees of Full Annular E? Combustor 


Viewing from casing toward hub Viewing from combustor exit 
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¢ One without a vane, two others with vanes set at different relative 
positions in relation to the fuel nozzle, “clocking”. 
¢ The difference of clocking between Case 1 and Case 2 is two degrees. , 


Tetrahedral mesh is generated by Cubit, and the total mesh count is 
approximately 50 million elements for all cases (AMR is off). 


Very uniform and fine mesh (~0.5 [mm]) inside the combustor and vane. 


Boundary Condition 


Aame Index Gas _ Mass flow rate {kg/s] 
Inflow Al Air 0.26 
Main dome Fl Fuel 0.00182* 
Ven! Pilot dome F2 ~~ Fuel 0.00182* 
eu : ———Lo a: Diffuser Bleed Sl Air - 0.018 
a (eS Pilot splash plate cooling A2 Air 0.0104 
ei 22 Pilot Dome A4 AS Outer liner cooling 1 A3 Air 0.0053 
mare irre d Outer liner cooling 2 A4 Air 0.0053 
Outer liner trim cooling A5 Air 0.0018 
Outer liner cooling 3 A6 Air 0.0024 
Main splash plate cooling | A7 Air 0.0116 
Inner liner cooling 1 A8& Air 0.0096 
Inner liner cooling 2 AQ Air 0.0056 
he gm Shingled tine Inner liner trim cooling AlO- Air 0.0018 
Construction Outer liner cooling 3 All Air 0.0024 
Counterrotating Centerbody outer cooling | Al2 = Air 0.0018 
=—* Centerbody mid cooling Al3 Air 0.0024 
Centerbody Inner cooling | Al4 _— Air 0.0024 
P3 [atm] T3(K] W3 [kg/s] Wrora [ke/s]f/a  Wepitor/Wrorar  Tyuet [KI] 
SLTO BY 4 720 0.26 0.00364 0.014 0.5 520 


Taken into consideration is the simulated sea level takeoff condition (SLTO), 
which is the most severe condition during the engine operation cycle. 


Cooling air is treated as source/sink terms on the surface 


Numerical Setup 


- LES-AUSM*-up 


¢ Liquid droplets (C12H21) are stochastically 
injected from the main and pilot domes with 
70° cone angle (hollow cone) 
¢ Finite-rate chemistry (2step-mechanism!")): 
KERO + 17.250, » 12CO,+ 10.5 H,O 
CO+050, » CO, 
ker = Anfi(pye or [KERO]™*° [02 }"°2", 
ky 2= A2f2()e! Eo2/RT) [co}"© [O2]"°22 , 
¢ Chemical integration: KIVA scheme. 
¢ Turbulence mode: non-linear k-e model with 
the wall function. 
About CPU 
¢ Used 1080 processors of Pleiades at NASA 
Advanced Supercomputing facility. 


¢ ~3 weeks to get statistics of unsteady 


calculations. 
[1] Franzelli, et a/., Combustion and Flame, Vol. 157, (2010) 10 


Numerical Results 


Red:100 m/s, Blue: -50m/s Red:100 m/s, Blue: -100m/s Red:100 m/s, Blue: -100m/s 
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Pressure Mach Number 


ill, 


¢ Captured eddies of many different length 
scales, central recirculation zone (CRZ) and 
precessing vortex core (PVC). 


Time-Averaged Flow Fields 


Case 2 


Case 1 


Case 0 
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Blocking Effect: Mach Number Sy 
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0 
¢ Due to the blocking effect, there are up to 20% difference in Mach 
number throughout the combustor. 
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Unsteady Temperature Fields (Case2) S 
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iiaeaain Temperature Fields @ Exit (P40) S 


P40: X=0.32[m] 
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Case 2 ms * timeimsy ; 
Example of Time-history T profile 


¢ Temperature profile is not uniform in the circumferential & radial directions. 
¢ Hot products are locally present, and the cooling air stays close to both the hub 
and the case. — 16 


Time-Averaged Temperature Profiles Ss : Ss 


Case 0 Case 1 


Case 2 
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Case 0 shows slightly lower temperature than that 
of Case1 and Case 2. 


¢ Not much effects from clocking 
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Time-Averaged Temperature Profiles, con S 
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Again, not much difference seen between Case 1 and Case 2 
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Due to the blocking effect, Case 0 shows lower temperature, especially 
in the main, but Case 1 and Case 2 are similar to each other. 
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Blocking Effect: Temperature Sy 


Clocking Effect on Temperature Fields 
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At P40, there are some noticeable 
differences between Case 1 and Case 2. 
For Case 1, there are five distinct hot 
regions, but only three for Case 2. 

Also, the two hot regions in the top parts 
of Case 1 are slightly hotter than the 
ones in Case 2 (more mixing?) 


Unsteady Temperature Fields @ Vane e& 
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¢ Adiabatic condition is used and any cooling holes on the surface are not included. 
¢ Locally hot (1600 [K]) and cool (900-1000 [K]) spots co-exist, which can cause the 


severe thermal stress on the vane, especially on the pressure side. a 


Time-Averaged Temperature Profiles @ Vane Sy 


Case 2 


¢ For Case 1 (top), hot streaks (1600 [K]) located at the far-right vane 
(V1), and these spots elongate on the pressure side. 

¢ For Case 2 (bot.), only one hot streak on V1. Also, the temperature of 
these hot streaks (1450 [K]) is not as high as the ones in Case 1. 


Histories of Temperature Profiles @ Vane 
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Time [ms] 
¢ Case 2 shows stronger temperature variation, with lower mean temperatures. 


(Standard deviation of temperature of v1: Case1 176 [k], Case2 200 [K]) 


Conclusions and Future Work 


In this paper, we have investigated the E° combustor using the 
OpenNCC developed at NASA Glenn Research Center. 


The three different geometries, without the vane and with the vane of 
two different clockings, are considered. 


Our numerical results show: 


LJ Presence of the vane at the combustor exit increases the pressure inside the 
combustor, resulting in weaker swirling flows and dilution airflows. This blocking 
effect causes up to 20% difference in Mach number. 

LI The effect of clocking is not significant in the flow fields inside the combustor and 
temperature fields in the primary zone. However, it indeed affects the distribution 
of the hot-streaks at the first stage vane surface as well as P40. 

LI In addition, one case shows stronger temperature variation with time, but cooler. 
Also, our results show that the temperature field on the pressure side is much 
higher than on the suction side. 


Future works will focus on the comparison of turbulence chemistry 
interaction models and adding the rotor. 
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Thank you! 


Questions? 
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